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Abstract 

In this article we analyze the geodesies of test particles and light in the five dimensional (charged) 
doubly spinning black ring spacetime. Apparently it is not possible to separate the Hamilton-Jacobi- 
equation for (charged) doubly spinning black rings in general, so we concentrate on special cases: 
null geodesies in the ergosphere and geodesies on the two rotational axes of the (charged) doubly 
spinning black ring. We present analytical solutions to the geodesic equations for these special cases. 
Using effective potential techniques we study the motion of test particles and light and discuss the 
corresponding orbits. 

1 Introduction 

A promising candidate to solve the problem of quantum gravity is the string theory. For its internal 
consistency more than four dimensions are required. This led to growing interest in higher dimensional 
solutions and especially higher dimensional black holes (see e.g. [I]). The higher dimensional analogon of 
the stationary axisymmetric Kerr black hole was found by Myers and Perry They also conjectured 
the existence of higher dimensional black holes with a non-spherical event horizon. Emparan and Reall [3] 
found such a solution with a topology of S 1 x S 2 , the rotating black ring. Since there is a certain range 
of values for the angular momentum where five dimensional black holes as well as two branches of black 
rings exist, the uniqueness theorem is no longer valid in higher dimensions. 

The first black ring solution found possesses a single angular momentum. A black ring with two angular 
momenta rotating in two independent planes was constructed by Pomeransky and Sen'kov [3]. 

In 2003 Elvang presented a charged singly spinning black ring solution [5] . The doubly spinning version 
of the charged black ring was found by Hoskisson [6] . Gal'tsov and Scherbluk [7] constructed a three-charge 
doubly spinning black ring, but unless the Dirac-Misner string is removed by setting two of the charges to 
zero, this solution is unbalanced. Various further black ring solutions are known (see e.g. j8] 19] ITU] 111]). 



The geodesic equations for test particles and light in Myers-Perry black hole spacetimes are separable 
[1211131 HPi] . Apparently this is not possible in black ring spacetimes, in general. However one can separate 
the equations of motion in special cases: geodesies on the two rotational axes (which are actually planes) 
of the charged doubly spinning black ring and zero energy null geodesies (which can only exist in the 
ergosphere) [15 1 IT6 ] . 

The geodesic motion in the singly spinning black ring spacetime was studied by Hoskisson [15] . He 
discussed the separability of the Hamilton- Jacobi equation and studied numerically the motion in the 
equatorial plane and the rotational axis of the black ring. Elvang, Emparan and Virmani [17] analyzed 
some aspects of null geodesies in the equatorial plane. Also zero energy null geodesies in the singly 
spinning dipole black ring spacetime were studied by Armas [18] . Igata, Ishihara and Takamori [19] found 
numerically stable bound orbits on and near the rotational axis. 

Analytical solutions of the equations of motion for zero energy null geodesies, geodesies in the equatorial 
plane and on the rotational axis of the singly spinning black ring were given in [20] , where also the possible 
orbits were discussed. 

Geodesies of the doubly spinning black ring were studied by Durkee [16] , who separated the equations 
of motion for zero energy null geodesies and on the two rotational axes and analyzed the effective potential. 



But so far neither numerical nor analytical solutions of the equations of motion have been given in the 
(charged) doubly spinning black ring spacetime. 

Hagihara [21) was the first to solve the equations of motion in a black hole spacetime analytically. 
He found the solution of the equations of motion in the Schwarzschild spacetime in terms of the elliptic 
Weierstrafi p function. 

When the cosmological constant is added to the Schwarzschild metric, the geodesic equations are 
of hyperelliptic type. The analytical solution of the geodesic equation in Schwarzschild- (anti) dc Sitter 
spacetimes was found in [22]. Also in higher dimensional Myers- Perry spacetime [23] and in higher 
dimensional Schwarzschild, Schwarzschild-(anti)de Sitter, Reissner-Nordstrom and Reissner-Nordstrom 
-(anti) de Sitter spactime [23] the equations of motion were solved analytically. 

Mathematically speaking the integration of the equations of motion can be reduced to the solution of 
the Jacobi inversion problem. The solution can be found by restricting the problem to the theta-divisor, 
the set of zeros of the theta function. This method was developed in 2003 by Enolski, Pronine and Richter 
[25[ to solve the problem of the double pendulum. 

In this article analytic solutions of the geodesic equations in the charged doubly spinning black ring 
spacetime are presented (for special cases). In the case E — m — the equations of motion are of 
elliptic type, however on the two rotational axes of the doubly spinning black ring, the equations are of 
hyperelliptic type. 



2 (Charged) Doubly Spinning Black Ring Spacetime 



The metric of an uncharged doubly spinning black ring can be written in the form [16] 
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The metric is given in toroidal coordinates (see [I} where — 1 < £ < 1, — oo < y < — 1 and — oo < t < oo. 
<f> and tp are 27r-periodic. The metric functions are 
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The parameters R, X and v describe the shape, mass and angular momenta of the ring and lie in the range 
< v < 1 and respectively 2^/i> <A<1 + za lf^ = the metric reduces to the singly spinning black 
ring. 

A ring-like curvature singularity is located at y = —oo. The metric has a coordinate singularity at 
G(y) — 0, so there are two horizons 
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The ergosphere is determined by H (y, x) — 0. In the singly spinning case v = the ergosphere has 
S 1 x S 2 topology like the horizon, but in the doubly spinning case v ^ the shape of the ergosphere 
becomes more complicated. Depending on the parameters v and A the shape of the ergosphere varies 
from ring- like S 1 x S 2 topology to spherical 5* 3 U S 3 topology (see [TB] for a detailed description) . There 
also exists an ergosphere- free region behind the inner horizon (y < yh-), which as not mentioned in [16] . 
An example for a ring-like ergosphere is shown in figure 2(a) and 2(b) and an example for a spherical 
ergosphere is depicted in figure 2(c) and |2(d)| . 



Figure 1: Toroidal coordinates (or ring coordinates) on a cross section at constant angles 4> and ip. Solid circles correspond to 
y — const, and dashed circles corrcdspond to x — const.. 1261 



The doubly spinning black ring has two angular momenta and thus the rotation is described by 

Q = Q^dip + Sl<f,d<f) , (5) 

where 
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and 
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One can split the polynomials A(x,y) and L(x,y) into x- and y-parts which will be useful later: 
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In [7] a three-charged version of the doubly spinning black ring was constructed. Here we use the 
balanced version where two of the charges are set to zero to remove the Dirac-Misner string. Compared 
to [7J the metric of the charged doubly spinning black ring is written in a slightly different way so that 
the metric functions G(x), H(x,y), A(x,y), L{x,y) stay the same as in the uncharged case from [TB] . 

The metric is 
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(a) v — 0.1 and A — 0.7: ring-like crgosphcrc in a-b coordi- 
nates 



(b) v = 0.1 and A — 0.7: ring-like ergo- 
sphere in X3-X4 coordinates 
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(c) v = 0.2 and A 
dinatcs 



0.9: spherical crgosphcrc in a-b coor- 



(d) v — 0.2 and A — 0.9: spherical crgo- 
sphcrc in coordinates 



Figure 2: Different shapes of the crgosphcrc plotted in a-b coordinates and X3-X4 coordinates. The a-b coordinates describe a plane 
of constant angles and ip. The £3-2:4 coordinates describe the equatorial plane of the ring where x — ±1. A more 
detailed explanation of the different coordinates and planes will be given in section [6] The crgosphcrc is depicted as a 
light red area with a red dotted border. There is no crgosphcrc in the white region. The black and grey dashed circles arc 
the inner and outer horizons of the black ring. 



The function D(x 1 y) is defined as 
D(x,y) = c 2 
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The charge is represented by the parameters 

c = cosh(a) and 
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The horizons have the same coordinates as in the uncharged case (see equation Q) and the ergosphere 
is still determined by H (y, x) = 0. If c = 1 and s — then D(x, y) — 1 and one obtains the metric of the 
uncharged black ring. 



The non- vanishing components of the inverse metric of the charged doubly spinnig black ring are 



9 H 


= -D(x,y) 2 




= cD(x,y)~ 




= cD{x,y)~ 


g H 


^Dix^y 1 




= -D(x,y)- 




= -D(x,y)- 


g xx 


= D(x,y)- 1 


gVV 


= -D{x,y)- 



/3 H(x,y) 2 1/3 {x-yf 

H(y,x) +C U{X ' V] R?H{x,y) 



R?H{x,y) G(x)G(y) 
1/3 ( X ~V) 2 ^<f,L{x,y) + n^A(y,x) 
R 2 H(x,y) G(x)G{y) 

1/3 ( x -y) 2 A{x,y) 



n\A(y,x) - 2n^L(x,y) - ^%A{x,y) 
G(x)G(y) 



R 2 H(x,y) G(x)G(y) 

/ 3 ( x -v) 2 My,*) 

R?H{x,y) G(x)G(y) 

/3 ( x ~yf L ( x ,y) 

R 2 H(x,y) G(x)G(y) 
1/3 (x-y) 2 (l-z/) 2 



(13) 



R 2 H(x,y) 
x — y) 2 (l — v^ 2 
R 2 H(x,y) 



G(x) 



-G(y) 



The metric of the charged doubly spinning black ring and its Hamiltonian — ^g ah PaPb do not depend 
on the coordinates t, <fi and "0, so we have three conserved momenta p a = g a bX b with the associated killing 
vector fields d/dt, d/d(j> and d/dip. A dot denotes the derivative with respect to an afhne parameter r. 
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E is the energy, $ and ^ are the angular momenta in <fi- an ^-direction. The conjugate momenta in x- 
and y-direction are: 
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To obtain the equations of motion for an uncharged particle in the charged doubly spinning black ring 
spacctime we need the Hamilton- Jacobi equation: 



dS_ 



0. 



(19) 



We already have three constants of motion (E, <& and ^/) and the mass shell condition g ab p a Pb = —m 2 
gives us a fourth , so we can make the ansatz 
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Inserting this ansatz into (|T9| gives 
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The Hamilton- Jacobi equation does not seem to be separable in general. However, it is possible to separate 
the equation in the special case E = m = 0. These zero energy null geodesies are only realisable in the 
ergoregion. 

We can also obtain equations of motion for geodesies on the <p- and ip-axis by setting y = — 1 (V>-axis) 
or x = ±1 ((/>-axis). The plane y = — 1 which is called the ip-axis corresponds to one of the two rotational 
axes of the doubly spinning black ring. The plane x — ±1 which is called the </>-axis, is the equatorial 
plane of the black ring and also the second rotational axis. 

In the next chapters we will study these special cases and solve the corresponding equations of motion 
analytically. 



3 Null geodesies in the Ergosphere 

For E = m = it is possible to separate the Hamilton- Jacobi equation: 
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With a separation constant k, the equation (|2"2"|) splits into two: 
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where 



U(x) = -[-/3(x)$ 2 - 26(x)<S>f + a(x)f 2 + kG(x)} (25) 
V(y) = -[a(y)^ 2 -2S(y)m-l3(y)^ 2 + kG(y)}. (26) 

Using p a = -j^ and (|T4| -(|18 p the separated Hamilton- Jacobi equation gives the equations of motion: 

dx 

dy 
ck, 



= ^W) (27) 



-VW) (28) 
[A{x,y)$- L(x,y)V] 



d 7 G(x)G(y) 

cb/> [-L(x,y)$-A(y,x)V] 

d 7 G(x)G(y) 



(29) 
(30) 



where 
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X(x) = -(l-v) 2 U(x) (32) 
Y(y) = -(l-u) 2 V(y) (33) 

are polynomials of fourth order. We also introduced the Mino-time [27. d'y — D(x, y)^ 1 ^ 3 ^ dr. 

From (|2T|) - (I5T|) we see that the x-, y-, <fi- and i/'-equation of motion do not depend on the charge. 
Therefore the charge of the black ring has no effect on the possible orbits for null geodesies in the ergo- 
sphere. 

In the sections below we will solve the equations ([27l) - (|28j) analytically. The i-equation (|3Tj) cannot be 
solved in general because it cannot be separated into x- and y-parts. Only in the special case v = which 
corresponds to the singly spinning black ring, it is possible to solve the i-equation (see |20j ) 



3.1 Classification of geodesies 



(1 - is) 2 Vd7 
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Equation (j2"7| and (|28p can be written as 

' ' d ''' ' + U(x) = (34) 

(!-*)■ v*r,' = °' (35) 

t/(x) and ^(y) can be regarded as effective potentials (see [IB]). To get real solutions for the x- and 
y-equation the effective potentials have to be negative. The zeros of the effective potentials (because 
E = 0) and hence X and Y mark the turning points of the motion of light or a test particle (in this case 
we only have light since m = 0). Thus the number of zeros determines the type of the orbit. 

Since X and Y are polynomials of fourth order it is difficult to calculate the zeros analytically. A good 
way to determine the number of zeros are parametric diagrams. 

For the x- motion we use the equations U(x) = and d = to construct a parametric ^-^-diagram. 
The diagram for the y-motion can be found analogously. If both $->I'-diagrams are shown in the same 
plot they form regions for the different numbers of zeros. Here one has to keep in mind that — 1 < x < 1 
and — oo < y < — 1 to find the allowed regions. 

Figure [3] shows the parametric ^-^-diagram for the x- motion (grey lines) and y-motion (black lines). 
Inside the grey and black coloured region different types of orbits are possible. Here U (x) has zeros in the 
range — 1 < x < 1 and V(y) has zeros in the range — oo < y < — 1. Outside the coloured regions U(x) has 
no zeros or is positive for — 1 < x < 1 and therefore no motion possible. 

If |$| > and \P lies in the coloured region the effective potential U(x) has always two zeros between 
— 1 and 1. If $ = the effective potential has the zeros —1 and 1 because 

U(±l) = (l-i/) 3 (l + i/±A) a $ 2 . (36) 

Also a third zero between —1 and 1 is possible in the case $ = 0. 

From (|36p it is clear that E/(±l) is always positive for |$| > 0, which means that there is a potential 
barrier preventing the photons from reaching the "equatorial plane" at x = ±1. One can also conclude 
that U{x) is negative between its two zeros and the x-motion takes place between these two values. 

Inside the grey region of the parametric diagram V(y) has two zeros. For small v and A a new region 



appears (the black region in figure 3(a) ), here V(y) has only one zero. 
At y = — 1 we have 

V{-l) = (l-v) 2 (l + v -A) 2 * 2 . (37) 

For ^ = one of the zeros lies at y = —1, in this case V(y) has up to three zeros. If |^| > a potential 
barrier prevents photons from reaching y = — 1 since V(— 1) > for l^l > 0. 

The influence of the separation constant k can be seen in figure 3(c) and |3(d)| If k becomes smaller 



the coloured region, where orbits are possible, of the 'f-^-diagram shrinks. For k = the coloured region 
vanishes, now only solutions with $ = \I r = are possible, and for k < no orbits are possible. 



Possible orbits for null geodesies in the ergosphere of a (charged) doubly spinning black ring are: 

1. Terminating Orbits (TO): 

The photon approaches the black ring, crosses both horizons and falls into the singularity. 

2. Bound Orbits (BO): 

The photon circles the black ring on a periodic bound orbit. In the ergosphere these orbits can only 
exist on the axis of ^-direction [y = —1), moreover a special condition for the separation constant 
has to be fulfilled. 

3. Many-World Bound Orbits (MBO): 

The photon circles the black ring on a periodic bound orbit, but crosses both horizons several times 
on its flight. Everytime both horizons are traversed twice, the photon emerges into another universe. 

From the effective potential U (x) we see that the x-motion is always bound. Either the motion takes 
place between two zeros or x is constant at —1 (this is possible if U(x) has three zeros in the special case 
$ = 0). To determine the type of the orbit we have to consider V(y), too. Depending on the number of 
zeros of V(y) and accordingly Y(y) three different orbit types are possible. Table [TJ shows an overview of 
all possible orbit types, some examples of the effective potential can be seen in figure HI [5] and M ■ 



(a) k = 1, 1/ = 0.1 and A = 0.7 (b) k = 1, v = 0.1 and A = 0.9 

For small v and A there are regions (black) where V(y) For larger A, V(y) has always two zeros, 
has only one zero. 




(c) k = 0.1, v = 0.1 and A = 0.9 (d) k = 0, v = 0.1 and A = 0.9 

If the separation constant k becomes smaller, the al- If k = the allowed region vanishes, 
lowed region also becomes smaller. 



Figure 3: Parametric 'I' -<I>- diagrams of the charged doubly spinning black ring for E — m — 0. In the allowed (coloured in grey and 
black) region U(x) has two zeros. V(y) has two zeros in the grey and one in the blue region. In the special cases <3> — 
and/or * — 0, three zeros of U(x) and V(y) are possible (sec text). 



Type A: 

Y(y) has one zero in the range — oo < y < — 1. X(x) has two positive zeros between x = — 1 and 
x = 1. Only TOs are possible. 

Type B: 

Y(y) has two zeros. Since the orbit crosses both horizons, MBOs are possible. 
Type C: 

This orbit type is only possible for ^ = 0. Y(y) has three zeros. Between two of the zeros a MBO 
can be found, the third zero is located at y = — 1 and belongs to a BO. X(x) has a positive and a 
negative zero. If additionally $ = the turning points of the MBO lie directly on the horizons since 
then V(y) = -cG(y). 

The BO lies on the axis of -0-rotation (y = —1). Since the shape of the ergosphere varies from 
ring-like to spherical, sometimes the axis y = —1 is not in the ergosphere and sometimes parts of 
the axis are inside the ergosphere. If the axis y = — 1 is not in the ergosphere, no such bound orbit 
should be possible. For this reason the separation constant k has to fulfill a specific condition if 

^..^ "'^"V"^' 2 ^" - (38) 

1 — A + V 

If k = fc ax i s then every BO at y = — 1 which would be lying outside the ergosphere is forbidden by 
X{x) (X will be positive and will have no zeros). A BO at y = — 1 in the ergosphere will still be 
allowed and its turning points will be located at the border of the ergosphere. 

fcaxis can be found if the x-equation for E — m = (|27[) is compared with the x-equation on the 
V'-axis (see section 0]). Both equations have to be the same for y — — 1 and E = m = 0. 
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Table 1: Types of orbits of light in the (charged) doubly spinning black ring spacetime for E — rn — 0. The thick lines represent 
the range of y and the turning points arc shown by thick dots. The horizons arc indicated by a vertical double line. The 
single vertical line at the loft end is the singularity. The coordinate y ranges from — oo to — 1. Since the BO of type C is 
located at y — — 1, the separation constant has to fulfill a certain condition (equation J38t V 



3.2 Solution of the x-equation 

Equation (|2"7|) can be written as 
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where X is a polynomial of fourth order with the coefficients 
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(a) Potential U{x) (b) Potential V(y) 



Figure 4: k = 1, v = 0.1, A = 0.7, * = -0.2 and * = 3.6 

Effective potentials U(x) (red) and V(y) (blue). The vertical black dashed lines indicate the position of the horizons. The 
horizontal green dashed line represents the energy and the green dot shows the position of the turning point. In the right 
picture an example of an orbit of type A can be seen. 




(a) Potential U{x) (b) Potential V(y) 



Figure 5: k = 1, u = 0.1, A = 0.7, $ = 0.4 and * = 0.8 

Effective potentials U(x) (red) and V{y) (blue). The vertical black dashed lines indicate the position of the horizons. The 
horizontal green dashed line represents the energy and the green dots show the position of the turning points. In the right 
picture an example of an orbit of type B can be seen. 



(a) Potential U(x) 



(b) Potential V{y) 



Figure 6: k = 1, v = 0.1, A = 0.7, * = 0.4 and * = 

Effective potentials U(x) (red) and V(y) (blue). The vertical black dashed lines indicate the position of the horizons. The 
horizontal green dashed line represents the energy and the green dots show the position of the turning points. In the right 
picture an example of an orbit of type C can be seen. The bound orbit at y — —1 is not allowed for the chosen parameters 
v and A, since U(x) has no zeros if one chooses the right value for the separation constant (equation J38l >). 

The substitution x = ±— + xo, where xq is a zero of X(x), reduces the polynomial X(x) to third order: 

b x ,3U 3 + b X:2 u 2 + b x .iu + b xfi . (45) 



i \ 2 

au 



6r) 



The coefficients are 



b x>3 = ±(Aa xA xl + 2>a x , 3 x% + 2a x ^x + a x ,i) (46) 

b x ,2 = 6a x ,4 x o + 3a X:3 x + a x>2 (47) 

b x ,i = ±(Aa xA x + a x , 3 ) (48) 

b x .o = a xA . (49) 

A further substitution u = -^-^ (^4v — transforms that into the standard Weierstrafi form 



3 
2 

d7 



^) =4v 3 -g x , 2 v 2 -g x , 3 := P x , 3 (v) (50) 



where 

g x , 2 = — -^- and g x , 3 = — - — . (51) 

Equation (l50|) is of elliptic type and is solved by the Weierstrafi elliptic function [28] 

v(l) = p(7-7in> 9x,2,g x ,3), (52) 
where 7/ = 7 in + P° —7= =5 dv' and u x>in = ^-(x in - xq)' 1 + -f#. 

°'- m • v / 4u' J -g x ,2«' 2 -Sx,3 * lz 

Then the solution of (|27|) takes the form 

b • 3 

z(7) = ± — ; —+*o- (53) 

4p(7-7i„, 9x,2, 9x,3)- ~t 

3.3 Solution of the y-equation 

Equation (|28|) can be written as 

2 



d"; 



Y(y) = a yA y 4 + a y . 3 y 3 + a^ 2 y 2 + a y ,iy + a yfi , (54) 



where Y is a polynomial of fourth order with the coefficients 



a yA — 


-(1- 


i;) 2 [*V(l - A 2 + u(l - v)) - 2V<P\v + <5> 2 v{l - A 2 ) - 


- cv\ 


(55) 


a V,3 = 


-(1- 


v) 2 [(i$ 2 - $ 2 )A^ 2 (3 - 2v) - Ac] 




(56) 


a y , 2 = 


-(1- 


v) 2 [^> 2 v(2\ 2 + v{\ - v)) + 2*$(A 2 i^ + A) + $ 2 (2A 2 ^ 4 








+c{v - 


-1)] 




(57) 


%,i = 


-(1- 


v ) 2 [-m 2 \(l + (1 - vf) + 2*$AV^(1 - v 2 ) + <5> 2 \v(2 


-2,v) 






+cA] 






(58) 


%,Q = 


-(1- 


i/) 2 [-* 2 (l + A 2 ) - 2*$A 2 V^ - $ 2 ^(A 2 - v 2 + v + 1) 


+ c]. 


(59) 



The problem can be solved analogously to the x-equation. The polynomial Y(y) is reduced to third order 
by the substitution y — ±i + y ("Ho 1S a zero of Y(y)): 



du x " 



where 



^ ) = ty^U 3 + by^U 2 + by^U + by , (60) 



b y ^ = ±(Aa yA yl + "ScLy^yl + 2a y , 2 yo + a y ,i) (61) 

b y ,2 = 6a yA y 2 + 3a V:3 y + a ya (62) 

b yA = ±(4a a! 4j/o + a y ,3) (63) 

b y ,a = a V A • (64) 



Next the polynomial is transformed into the Weierstrafi form by substituting u = -r— ( 4u — 



3 



dv N 2 



d') 



4w 3 - 5y,2« 2 - 5 „, 3 := P y ,s( v ) > (65) 



where 



& a,2 b y^ b v,3 QT i . by A b ya by3 b yfi b 2 yZ b y2 

9y < 2 ~l2-—T- ^ 48 16~~216- (66) 

Since ([63)1 is solved by the Weierstrafi p-function, the solution of (|2"5)) yields 

»(7) = ±— r-^ + (67) 

4p(7-T&»^,a,5y,3)- T 

where < = 7in - J^ -^—i— — dt/ and ^, in = ^(y in - yo)" 1 + %f ■ 
3.4 Solution of the 0-equation 

To solve (f2"9")) the x- and y-parts have to be separated. Inserting A(x,y) and L(x,y) from ([5]) yields 

d0 _ /?(x)$ + «(?/)$ + <5(y)* 

d^ = i? 2 G(x) + R 2 G(y) ' 

Using (27) and ((28]) we get 

/3(x)$ + da; a(y)5> + <5(y)\& dy 



i? 2 G(x) ^) i? 2 G(y) v ^)- 

So the ^-equation consist of a x-integral and a y-integral: 

rx p(x')$ + 5(x')V dx' f y a(y')<S> + %')* dy' 



(69) 



r 2 g(x') 4„ i? 2 G(y') v/yry) 

=: + (70) 



Let us first consider I x . Here we substitute x = ± — J b 3 + xo to transform the polynonial X(x) into the 

iu - js i- 

WeierstraB form P x 3: 



4= / — ^ 2 7^ " T J " 7==- (71) 

4«' 

This integral has four poles p\ , P2, P3 and p4 (zeros of the function G with respect to u) . We next apply 
a partial fractions decomposition upon (|T1|) : 



(72) 



ifj are constants which arise from the partial fractions decomposition and depend on the parameters of the 
metric and the test particle. Then we substitute u = p(v, g x ,2,g x ,3) with p'(v) = -\/4p 3 (w) — g x ,2p(v) ~ 9x,3- 
Equation (|72"|) now yields 



4 K, 



I x = ^o + E^V ) dv ( 73 ) 



With V x = V x (j) = 7 - 7i' n and v x,in = "x(7in)- 

The second integral ij, can be rewritten by canceling the factor (1 — y 2 ): 



I, 



u[-(l + A 2 ) - - 1/) + Ay' (2 - 3i/) - (1 - A 2 )?/' 2 ] • $ + A^[A - (1 - ^ 2 )y' - A^y' 2 ] ■ * dy' 

i? 2 (l + Ay' + ^' 2 ) 

(74) 

We substitute y = ± — S£ ^-=- + yo to transform the polynonial Y(y) into the Weierstrafi form P y 3 and 
4« — 

apply a partial fractions decomposition upon (f?4"|) where the constants Lj arise. Note that the two poles 
Pi and P2 of / y are also poles of I x . 



2 La \ du' 



Iy= Xq + V^ . (75) 

If we now substitute u = p( v > 9y,2, gy,3) with p'(v) = -\/4p 3 (w) — g Vi 2p(v) — g y .3, equation (|72t yields 

-f Py{v)-Pj 

with = u y (7) = 7 - 7;" and u y , in = %(7in)- 

After solving the elliptic integrals of the third kind in ([73| and (|76l) (see e.g. [29]), the final solution is 

0(7) = <^ IfoK - Ux,in) + V , r \ 2 <x{Vj){v x - V xM ) + log f J -\ - log ■ 

+ \ L Q (v y - v v ,- m ) + V -y-^— f 2Cy(u i )(u J/ - v yM ) + log — Hl^ - log ■ 
1 J=l S'yW) V ^K.in-Wj) 



; log — 

+ W j) 

V j) ^ Vy(Vy, in +Vj) ) 



(77) 



with pj — p(vj) and 

Px(v) = p{v,g x ,2,gx,3) , p y (v) = p(«, g y ,2, g y ,3), 

Cx(v) =C(v,g x ,2,gx,a), ( y (v) = C(v,g y ,2,g v ,a) , (78) 

<r x (v) = <t(v, g x ,2, 9x,3) , <r y (v) = cr(v,g y ,2,g y ,3) ■ 



3.5 Solution of the ?/>equation 

To solve (f5Uf the x- and y-parts have to be separated. Inserting A(x,y) and L(x,y) from ([SJ yields 

dip _ 5(x)$ - a(x)* 6(y)® + P{y)^ 
ch7 ~ R 2 G{x) R 2 G{y) ' 



(79) 



Using d2TJ) and J28J) we get 



or 



5{x)<S> - a(x)V dx 5{y)<5> + 

cv ' i? 2 G(x) ^fo) v 7 ^) 

^ ^ _ f x 5(x>)$ - a(x')^> dx' + dy' (gl) 



i? 2 G(.x') J yin R 2 G(y') ^YW) 

The -0-equation can be solved analogously to the </>-equation from the previous section. Here the solution 



is 

( 2 



^( 7 ) = < M («x " ^4n) + V "77^7 2Cx(«i)(«x " «x.i„) + log 3 \ - log 



~[ Px\ v j) V CT x (^ !in - Uj) Cr x (« aj i n + Uj) 

-v„i <„ <■,,,„ i + V — ^- ( 2,, ( r, 1 1 ,„ - ,„,„ i + log "v^v-^ - log . a ^ + ^ ^ 



Vy(Vy,in ~ Vj) Vy(Vy,in + Vj) 



(82) 



where v = v{^) = 7 - 7" , v in = f(7in)) Pj = p(^) and 



Cx(v) =^{v,g x ,2,9x,3), ( v (v) =C(v,g y ,2,g v ,a), (83) 
CxCu) =(r(v,g x ,2,9x,a), v v {v) = g y ,2,g y ^) ■ 

Mj and TVj are constants which arise from the partial fractions decomposition and depend on the para- 
meters of the metric and the test particle. 

4 Geodesies on the -0-axis 



The surface y — — 1 is the axis of ^-rotation of the doubly spinning black ring. Here the Hamilton- Jacobi 

as 
'hi 



equation depends on the coordinate x only. We set y = — 1, ^ = and p y = = in the Hamilton- Jacobi 



equation (pH) : 

. ra ,. C V3 (l ,_ 1) ^ £S + D -V3 (l ,_ 1) (£±|^{ G(l) ( 
($ + cft £) 2 [/3(x) i/[2 + - 1/) + A(2 - 3i/)] 



dx 



(1-i/) 2 LG(s) 1-A + v 

This can be rearranged to 

'5SA 2 _ DVs^, -1)^(3;, _i) / H{x, 1) 



(84) 



{x + l) 2 {l-v) 2 G{x) 

i^ + cQ^E) 2 ( f3{x) u[2 + v{\ - v) + A(2 - 3u)} 



{l-v) 2 G{x) \G(x) 1 - A 



D 2 ^(x,-l)^-^-E 2 - m 2 ) 
H(-l,x) ) 

=-.x s . 



(85) 



Then we have 



S = ^m 2 r - Et + $0 + / v A s < lr ■ 1 Mi ) 



Now we set the partial derivatives of S with respect to the constants m 2 , E and <f> to zero in order to 

(g±l) , 

D(x,-l) 1 / 3 R 2 H(x,-l)H(-l,x) 



obtain the equations of motion. With the Mino-time [37] = D , _ 1 ^i/3P2jj(2 -iigr-i x) ^ T *' le ec l ua tions 



of motion take the form 
d.r 



— = ^X{x) (87) 



d 7 (1 + v - \) 2 G(x) 

dt 2 D(x,-l)H 2 (x,-l) (x + l)H(x,-l)H 2 (-l,x) 

d^ - RE W+T) (l + v-\) 2 G(x) <* l *(* + <*l*E) (89) 

where 

X(x) = {1-v) 2 H{x 1 -1)H{-1,x){r 2 G{x) D{x,~l)H{ Xl -l)E 2 - D 1/3 (x,-l)H(-l,x) 



m 2 



{i-x + uY 



H(-l, x)$ + cRXy/v y/2{(l + vf - A 2 )(l - x 2 )E 



(90) 



and 



H(-l,x) = (I- X) 2 -v 2 + vx 2 (l- X 2 -v 2 +2Xv) 

H(x,-1) = 1 + A 2 - v 2 -2Xv{\- x 2 ) + 2A.t(1- v 2 ) + x 2 v{l- X 2 - v 2 ) 

D(x,-1) = 1 + ' _ 1} [2A(1 - u){x + 1)(1 + vx)\ 

RX^2{{l + v) 2 -X 2 ) ^ 

n * = — hf^x) — a-*)^- ( 9l ) 

Solving the equations of motion (|87p - (j8"9")l analytically is only possible if X(x) is a polynomial. This 
happens in two cases: 

1. D(x, —1) = 1 (which implies c = 1 and s — 0): uncharged doubly spinning black ring 

2. m = 0: charged doubly spinning black ring with photons 

In both cases X(x) is a polynomial of 10th order, so that the equations of motion are of hyperelliptic type. 

4.1 Classification of Geodesies 

From ([57)) we can read off the effective potential consisting of two parts U+ (x) and U- (x) : 

X = a(x)(E - U+)(E - U-) . (92) 
Since X(x) can be written as X(x) — a(x)E 2 + b(x)E + c(x) the effective potential takes the form 



TT , v ~b(x) ± y/b(x) 2 - Aa{x)c{x) 

U±{x) = v , 93) 

2a(x) 



where 



/ c 2 £l 2 H 2 (-l x)(x + I) 2 \ 
a(x) - (1 - v) 2 H(x, -1)H(-1, x) lR 2 G(x)D(x,-l)H(x, -1) - \ 1 _ x+ \ /)2 ~ 1 

b(x) = -2(1 - y)2$C n,g3 ( 1)(iC + 1) 2 (Q4) 



z(x) = -(1 - v) 2 H(x, -1)H(-1, x) (r 2 G{x)D 1/3 (x, -1)H(-1, x)m 2 + 



(x + iyH 2 {-i,x)<i> 

(1-A + ^) 2 



Figure [7] shows the effective potential for the motion on the ip axis. U+ is plotted in red (solid line) 
while U- is plotted in blue (dotted line). The grey area between the two parts of the potential is a 



forbidden zone where no motion is possible because X{x) becomes negative there. If $ — then X(x) 
has none or one zeros and U+ and U— are symmetric ([/+ = —U-). If |$| > there is a potential barrier 
which prevents the geodesies from reaching x = +1. X(x) can have up to three zeros. 

Possible orbits are Bound Orbits (BO), where light or test particles circle the black ring, and Escape 
Orbits (EO), where light or a test particle approaches the black ring, turns around at a certain point and 
escapes the gravitational field. 

There are five different types of orbits (see table which exist in the charged as well as in the 
uncharged black ring spacetime: 

• Type A: 

X(x) has no zero. EOs without a turning point exist. The orbit crosses the equatorial plane (x = +1) 
and reaches infinity (x = — 1 and y = —I). 

• Type B: 

X(x) has one zero. BOs with a turning point on each side of the ring exist, so that the orbit crosses 
the equatorial plane (x = +1). 

. Type C: 

X(x) has one zero. EOs with a turning point exist. 

• Type D: 

X(x) has two zeros. BOs which do not cross the equatorial plane exist. 

• Type E: 

X(x) has three zeros. BOs which do not cross the equatorial plane and EOs exist. 




(a) * = (b) $ = 1 (c) 4> = 8 

Examples of orbits of type A and B. Examples of orbits of type C and Example of an orbit of type E. 

D. If |$| > there is a potential The potential can have local cxtrcma 

barrier which prevents the geodesies which lead to three zeros of X(x) 
from reaching x — +1. X(x) can have 
one or two zeros. 



Figure 7: R = 1, m = 1, s = 1, v = 0.1 and A = 0.7 

Effective potentials U-^(x) (red solid line) and U-{x) (blue dotted line) on the i/;-axis (y — —1) of the black ring. The 
grey area is a forbidden zone, where no motion is possible. Green dashed lines represent energies and green points mark 
the turning points. 



4.2 Solution of the x-equation 

In the two cases D(x, —1) = 1 (i.e. s = 0) and m — equation (1571) can be written in the form 

/j \ 2 10 

i^-j =X{x)=Y^a i x i . (95) 

A separation of variables gives the hyperelliptic integral 

f x da;' . . 



type 


zeros range of x 


orbit 


A 


U 


' EO 


B 


1 • 


BO 


C 


1 


• ' EO 


D 


2 • 


• ' BO 


E 


3 * * 


• ' EO, BO 



Tabic 2: Types of orbits of light and particles in the charged doubly spinning black ring spacctimc for y — —1, = 0. The thick 
lines represent the range of the orbits. The turning points arc shown by thick dots. The coordinate x ranges from —1 to 
+ 1. An orbit of type B crosses the equatorial plane. 



Since X(x) is a polynomial of 10th order, the genus of the Riemann surface is g = 4. A canonical basis of 
holomorphic (du^ and meromorphic (dr,) differentials associated with the hyperelliptic curve w 2 — X(x) 
is given by (see [24] or [29] ) 



dui :— dx . 

w 



2g+l-i 



dn:= ^2 (k+l-i)\k+i- 



' Aw 



-dx 



(97) 
(98) 



with i = l,...,g and being the coefficients of the hyperelliptic curve written as 

w 2 = X 2g +2X 29+2 + X 2g+ ix 29+1 + X 2g x 29 + . . . + A . (99) 
Furthermore we introduce the holomorphic and meromorphic period matrices (2w,2u/) and {2rj 1 2r]'): 



= f dui, 


Mi ■■= 


f dm, 




Jbj 


- f dn, 


H 3 ■■= 


f dr 4 







(100) 



with i,j — 1, ...,g, where {aj, bi\i = 1, ■■■,(?} is the canonical basis of closed paths. 

The solution of equation (pM]) is extensively discussed in [25], [23], [23], [25] and is given in terms of 
the derivatives 



ai(u) = --^—a(u) 
aui 

aij{u) = - ——a(u) 
oui ouj 



(101) 



of the Kleinian sigma function a(u) = fc e -(i/2)« t ')" 1 «,9((2cj)- 1 m + K Xin ; r) : 



c(7) = _^e), 



where 



7e 



(102) 



(103) 



/ 7-Vin \ 
72 
73 

V 74 J 

The constant j' in = y ul + J x dui depends on y ln and x- ln only. 72, 73 and 74 are defined by the vanishing 
condition of the Kleinian sigma function: er(7 e ) = 0, 174(70) = an d cr 44(7e) ~ so that (2o;) _1 7 e is 
element of the theta divisor (the set of zeros of the theta function). 



an 



4.3 Solution of the 0-equation 

With {ST} equation ||8SJ) yields 

(x + l)i/(x,-l)iJ 2 (-l,x) 



d.r 



(1 + 1/ - A) 2 G(a;) 



This can be rewritten as 



Pq (x) dx 



($ + ctt^E) ^ = . (104) 

— Pawyxw' (105) 

where Pe(a^) is a polynomial of sixth order and Psix) is a polynomial of third order. So we have to solve 
the integral 

r x P 6 (x') dx' 



which has poles at pi t 2 



,W) v^o (1 ° 6) 

-a±va 2 -4 £ an j p 3 = \_ We apply a partial fractions decomposition upon (|106j) : 

3 7 



2// 



da;' 



(107) 



where ifj are constants which arise from the partial fractions decomposition. 

If D(x, — 1) = 1 (i.e. s = 0) or m = then X(a;) is a polynomial of tenth order and equation (|107[) 
consists of hyperelliptic integrals of the first and third kind. 



First we will solve the holomorphic integrals of the first kind f x , x dx' (in this case i = 0, 1, 2, 3). 

We introduce a variable u so that w — wo = f x , dx The inversion of this integral yields x(v) — — y\ 
(see section or [29 ), where 

/ w - w \ 

, /i(0) = / 2 (0) = / 3 (0)=0. (108) 



u = 2L 



/ w - w ^ 
/i(w - w ) 
/ 2 (w - w ) 
V h{v-v Q ) j 

The functions /i(w — Wo), f2(v — ^o) ancl /3( u — v o) can be found from the conditions a(u) — 0, 174(11) = 
and 044 (it) = 0. 

The vector 21^ identified with each branch point is defined as [21?] 



a* = 



du = liiJEi + 2uj e,;, i = l,. 



,10, 



with the vectors and € R 4 with entries equal | or 0. The matrix 

m = ' " x 

is called the characteristic of a branch point ej 
Equation (|107[) now reads 







••4 ] 




\ en 


• ■ £i4 / 



(109) 



(110) 



E 



Ki 



dx' 



~{ J Xln X> - Pi y/Xlx 7 ) 



+ K A (v - w ) + K 6 fx(v - w ) + K 6 f 2 (v - w ) + K 7 f 3 (v - w ) . (Ill) 



The remaining hyperelliptic integrals of the third kind can be solved with the help of the following 
equation (see |29|). 



W 



1 dx 

/ x — Z w 



rP T / r(z,W) 

=2 / du T (x,y) C / 

JP' I \J(e 2 ,0) 



e 2 ,0) 
(Z,W) 



du + K 00 ) - 2( V 'e' + V e) - \${Z, W) 



In 



111 



(112) 



P and P' are points on the hyperelliptic curve, Z is a pole, W = w(Z) and w 2 — X(x). The zeros of 
w 2 (:r) are the branch points of the curve w 2 . du is the vector of the holomorphic differentials of the first 
kind dui — ^-^-dx with i = 1, g. C an d a are Kleinian functions and is the vector of Riemann 
constants. 

The gth component (in this case genus g = 4) of the vector 3(Z, W) is 3 9 (Z, W) = and for 1 < j < g 
we have 

W 

3 j (z,w) = TTg 



9-3-1 



£ (-l) 9 -^+i^5 s - fc - J --x(e). 



(113) 



k=0 



The 5fc(e) are elementary symmetric functions of order fc built on g — 1 branch points e4, . . . , e 2g : Sq = 1, 
S 1 ! = e 4 + . . . + e 2g , etc. 



Finally the solution of the ^-equation ([88[1 is 



-i„ + ifi 



' du 



(Pi,W<) 



(e 2 ,0) 



du + J - 2(t 7 '£ / + - ~3(Pi, Wi) 



q(iy 2 (x)) ^(V^ 2 (^n)) 
n a(^(x)) n a(^(xi„)) 



+ Ka(v - vo) + K 5 fi(v - v ) + K 6 f 2 (v - Wo) + K 7 f 3 (v - v ) , 

(114) 



where = y^fe) and W^ 2 (x) = /* du ± f ( ( e P ;J l) du - K 



4. 4 Solution of the t-equation 

With (|HI|) equation ([89]) yields 



di = i? 2 £ 



,£>(£, -l)i? 2 (a;, -1) (a + l)H(x, -1)H 2 (-1, x) 



(x + 1) 



This can be rewritten as 



dt 



(1 + v - X) 2 G(x) 
P 4 {x) P 5 {x)\ dx 



cQ^,($ + dl<f,E) 



dx 



t-ti, 



x+1 P2(x)J ^/xjxj 
) has the zeros ; 

P A {x') , P s (a; / ) 



(115) 



(116) 



where Pd{x) are polynomials of order d. P 2 (x) has the zeros = A± ^ 2 — ■ So we have to solve the 
integral 



dx' 



l XiD \x' + l P 2 (x')J JX(tf) 
We apply a partial fractions decomposition where the constants Mj and Nj arise: 



t-*i, 



dx' 



(117) 



(118) 



If D(x, — 1) = 1 (i.e. s = 0) or m = then X(x) is a polynomial of tenth order and the equation (|118p can 
be integrated analytically. The hyperelliptic integrals of the first and third kind can be solved as shown 
in section 14.31 The solution of the ^-equation (|89D is 



t 



M 1 



du 



(-1,^X7=1)) 
(ea,0) 



1. 



du + K x J - 2(Ve' + ije) - -3(-l, 



q(^ 2 (x)) a(W^Q) 

+ U1 <t(W^)) n ^(^ 1 (^n)) 



+ M 2 (« - v ) + M 3 fi{v - v ) + Mif 2 (v - vo) + M 5 f 3 (v - v a ) 



2 



du 



(e 2 ,0) 



dit + K c 



2(Ve , + i 7 e)--3(Pi,Wi) 



q(^ 2 Qr)) a(W 2 Qr in )) 
" n a(Wi(x)) n (7(^1 (x in )) 



+ N 3 (v - v ) + Nifi(v - v ) + N 5 f 2 {v - v ) + N e f 3 {v - v ) 



+ th 



(119) 



where W t = y^fe) and W l > 2 {x) = J^du± f^ i} du - K 



5 Geodesies on the 0-axis 

The surface x = ±1 is the axis of <f> rotation. It can be seen as the equatorial plane of the black ring, 
which is divided into two parts. The first part x = +1 is the plane enclosed by the ring (or more precisely: 
enclosed by the singularity), which we will refer to as "inside" the ring. The second part x = — 1 describes 
the equatorial plane around the black ring (or more precisely: around the singularity), which we will refer 
to as "outside" the ring. 

If we set x = ±1, $ = and p x = ^| = in the Hamilton- Jacobi equation (j2"Tj) . it depends on the 
coordinate y only: 



o . m > - tf/ 3(±1 , rigEMtf - d-/3 (±1 , „ (*' - »W - "Y G(y) (w 



2 

{V + cCl^E) 2 \/3(y) i/[2 + i/(1-i/)=fA(2-3j/)]" 



This can be rearranged to 



G(y) 1 ± A + v 



(120) 



dSV D 1 / 3 (±l,y)R 2 H(±l,y) ( 2 „ /3 ^ i ,H(±l,y) 



m — D ' (±1 iA — — E 

ih,) {±l-y) 2 {l-v) 2 G{y) V { ,V> H{y,±lf 

(g + ca 4 ,E) 2 f p(y) _ I /[2 + I /(1- I /)tA(2-3i/)] \ 

{l-v) 2 G{y) \G{y) l±\ + u 1 " 



(121) 



S = \m 2 T -Et + Vip+ I yjY s dy . (122) 



Then we have 

.<? = , 
2 

Now we set the derivatives of S with respect to the constants m 2 , E and *f> to zero in order to obtain the 
equations of motion. With the Mino-time |27j drf = D ^ ±1 Jy/i f^ H y ) ^ T ^ ne equations of motion take 
the form 

dy 



where 



^ (±^-y)H{±l,y)H{y,±l) ^ 

^ = — ji+v±wm — (* +c ^ ( l24 ) 

df r)2lp D(±l,y)H 2 (±l,y) (±1 - y)H(±l, y)H(y, ±1) 

d^ = flS (±l- y )ff(y,±l) (l + ^iATOy) + ( 125 ) 



Y(y) = (1 - ^) 2 ^±1) {R 2 G(y){D 1/3 H(y, ±l)m 2 - M(±l, y)E 2 ] 
{±l-yf r 



(1± A + ^) 2 



(», ±1)* - ci?AV2((l + ^) 2 -A 2 )- 



1 + y 



(l-X + u) 



•(1 + A - v + yv(l - A - v) ± 2v{l - y))E 



(126) 



and 



H(±l,y) = (1 ± A) 2 - v 1 + y 2 v{\ - A 2 - v 2 qp 2A^) 

#(y,±l) = l + A 2 -^ 2 ±2Ai/(l- 2 / 2 )+2Ay(l-^ 2 )+y 2 ^(l-A 2 -!/ 2 ) 



,s 2 



2?(±l,y) = 1 + , [2A(1 - v){±\ -y)(lT vy)] 

a* = - RW2 ^ ( + Z - A2) ^-(i + x-» + Mi - a - ») ± Mi - V)) ■ d27) 

ii(y, ±1) l-A + fc' 



Solving the equations of motion (|123p - (|125[) analytically is only possible if Y(y) is a polynomial. This 
happens in two cases: 

1. D(±l,y) = 1 (which implies c = 1 and s = 0): uncharged doubly spinning black ring 

2. m = 0: charged doubly spinning black ring with photons 

In both cases Y(y) is a polynomial of 6th order, so that the equations of motion are of hyperelliptic type 
(genus g = 2). 

5.1 Classification of Geodesies 

From (|123p we can read off an effective potential consisting of the two parts V+(y) and V-(y): 

Y = a(y)(E~V+)(E~V-). (128) 
Since Y(y) can be written as Y(y) = a(y)E 2 + b(y)E + c(y) the effective potential takes the form 



2a{y) 



where 



I , c 2 n 2 H 2 (y,±l)(±l-y) 2 \ 
a(x) = (1 - ufH{±l, y)H(y 7 ±1) f -R 2 G(y)D(±l, y)H(±l, y) (1 ± A + I) 2 J 

6(x) = -2(1 - v) 2^HHy,±l)H(±l,yK±l-y)- (130) 
w v ; (1 ± A + j/) 2 

c(x) = (1 - (±1, y)ff (y, ±1) ^(^^(il, j,)F(„, ±l)m 2 - ^-^^0^. 

The two cases x = +1 (geodesies inside the ring) and x = — 1 (geodesies outside the ring) have to be 
discussed separately. 

5.1.1 Geodesies outside the ring 

Let us first take a look at the motion on the surface outside the black ring. Here we have x = —1. Figured] 
shows the effective potential V(y) for different values of the parameters. V+ and V- meet at the horizons. 
Mainly the angular momentum "J defines the shape of the effective potential, the charge parameter s has 
less influence. For "J = the potential is symmetric with respect to the y-axis (V(y) = 0) and Y(y) has 
one or two zeros. A potential barrier prevents the test particles and light from falling into the singularity 
If l^l > the potential is no longer symmetric and up to three zeros of Y(y) are possible. If \^\ is large 
enough a gap in the potential barrier appears which allows test particles and light to reach the singularity. 

Possible orbits outside the black ring are: 

1. Terminating Orbits (TO): 

Light or test particles approach the black ring, cross both horizons and fall into the singularity. 

2. Many-World Bound Orbits (MBO): 

Light or test particles circle the black ring on a periodic bound orbit, but cross both horizons several 
times on their flight. Everytime both horizons are traversed twice light or the test particles emerge 
into another universe. 

3. Escape Orbits (EO): 

Light or test particles approach the black ring, turn around at a certain point and escape the 
gravitational field. 

4. Two-World Escape Orbits (TEO): 

Light or test particles approach the black ring, cross both horizons turn around at a certain point, 
cross the horizons a second time and escape the gravitational field. Since both horizons are traversed 
twice light or the test particles emerge into another universe. 



(a) s = 0, 5 and * = 5 (b) s = 1 and * = (c) s = 1 and * = 10 

Examples of orbits of type A and Examples of orbits of type C and Examples of orbits of type D and 
B. Co Do 



Figure 8: R = 1, m = 1, u = 0.1 and A = 0.7 

Effective potentials V^(y) (red solid line) and V- (y) (blue dotted line) on the 0-axis outside the (charged) doubly spinning 
black ring. The grey area is a forbidden zone, where no motion is possible. The horizons arc marked by vertical dashed 
lines. Horizontal green dashed lines represent energies and green points mark the turning points. 

There are four different types of orbits (see table [3]) . 

• Type A: 

Y(y) has no zeros and only TOs exist. 

• Type B: 

Y(y) has one zero and only TEOs exist. 

• Type C: 

Y(y) has two zeros and only MBOs exist. In a special case the two turning points lie on the horizons. 

• Type D: 

Y{y) has three zeros. MBOs and EOs exist. In a special case the two turning points of the MBO 
lie on the horizons. 
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-• * MBO 
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3 ' 1 f 


• MBO, EO 



Table 3: Types of orbits of light and particles in the (charged) doubly spinning black ring spacctimc for x — —1 and $ — 0. The 
thick lines represent the range of y and the turning points arc shown by thick dots. The horizons arc indicated by a vertical 
double line. The single vertical line at the left end is the singularity. The coordinate y ranges from -co to —1. 



5.1.2 Geodesies inside the ring 

The effective potential for geodesies on the surface enclosed by the black ring {x = +1) is shown in figure 
[9l Again, if we have = the potential is symmetric and Y(y) has none or one zeros. In the case > 
a potential barrier appears which prevents test particles and light from reaching y = — 1. In contrast to 
the motion of light and test particles outside the black ring, the charge parameter s has now significant 
influence on the effective potential and the possible orbits. If s > up to four zeros of Y(y) are possible, 
if s = a maximum of three zeros is possible. In the case of four zeros, bound orbits of test particles (not 
light) behind the inner horizon are possible. 




(a) v = 0, 1, A = 0.7, m = 1 and * = (b) f = 0, 2, A = 0.9, m = and * = 5 (c) i> = 0, 2, A = 0.9, m = 1 and * = 
Examples of orbits of type A,C, D Examples of orbits of type B,E and 20 

and Do- Eo Examples of orbits of type F. 

Figure 9: R - 1 and s = 0, 1 

Effective potentials (y) (red solid line) and V_ (y) (blue dotted line) on the 0-axis inside the (charged) doubly spinning 
black ring (x — +1). The grey area is a forbidden zone, where no motion is possible. The horizons are marked by vertical 
dashed lines. Horizontal green dashed lines represent energies and green points mark the turning points. 



Possible orbits are Terminating Orbits (TO), Many- World Bound Orbits (MBO) and Bound Orbits 
(BO). The bound orbits, which are only possible in the charged case, exist behind the inner horizon of 
the charged doubly spinning black ring. 

There are six different types of orbits (see table 2]). 

• Type A: 

Y(y) has no zeros and only TOs exist. Test particles or light come from infinity (y — —1) and fall 
into the singularity (y = -co). 

• Type B: 

Y (y) has one zero and only TOs exist. The orbits starts at a certain point and ends in the singularity. 

• Type C: 

Y{y) has one zero and only MBOs exist. Light or a test particle crosses the center of the ring 
(x = +l,y = -l). 

• Type D: 

Y(y) has two zeros and only MBOs exist. In a special case the two turning points lie on the horizons. 

• Type E: 

Y{y) has three zeros. MBOs and TOs exist. In a special case the two turning points of the MBO 
lie on the horizons. 

• Typ F: 

Y{y) has four zeros. MBOs and BOs exist. The BOs lie behind the inner horizon and are only 
possible if the doubly spinning black ring is charged and m/0. 



5.2 Solution of the y-equation 

If Z?(±l, y) = 1 (i.e. s = 0) or m = then Y(y) is a polynomial of sixth order and equation (|123[) can be 
solved analogously to ([57)1 . A separation of variables yields 

f v dy' 

7-7m = -/ r—— ■ (131) 



The solution is 

( } = _^2e) . (132) 

02(76) 



where 

- ( - 

72 



7eH )=<Poo- C' : W) 
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Tabic 4: Types of orbits of light and particles in the (charged) doubly spinning black ring spacctimc for x — 4-1 and <& — 0. The 
thick lines represent the range of y and the turning points arc shown by thick dots. The horizons arc indicated by a vertical 
double line. The single vertical line at the left end is the singularity. The coordinate y ranges from -co to —1. 



The constant j' in = y ln — J°° , dy depends on 7; n and y- ln only. 72 is denned by the vanishing condition 

of the Kleinian sigma function: er(7 e ) = so that (2cj) _1 7q is an element of the theta divisor (the set of 
zeros of the theta function). 



5.3 Solution of the ^-equation 

With ([123| equation ([T24"| yields 

(±l-y)H(±l,y)H(y,±l) 



dtp 



This can be rewritten as 



(1 + v ± A) 2 G(y) 

Pi{y) dy 



dy 



Vny) 



dtp = 



P 3 (y) 0%) 



(134) 



(135) 



where Pi{y) is a polynomial of sixth order and Pz{y) is a polynomial of third order. So we have to solve 
the integral 

rv P A (y>) dy' 



(136) 



which has poles at pi^ = A± ^ 2 — = yn± and p 3 = —1 (if x = +1) or p 3 — 1 (if x = —1). We apply a 
partial fractions decomposition upon (I136p : 



-0 - Ipi, 



Ki 



\i=i 



y - Pi 



K 4 + K 5 y' 



dy' 



(137) 



where K{ are constants which arise from the partial fractions decomposition. 

If D(±l,y) — 1 (i.e. s — 0) or m — then Y(y) is a polynomial of sixth order and equation f|l 37[) 
consists of hyperelliptic integrals of the first and third kind. The solution can be found in the same way 
as in section 13.41 or 14.31 



i> = V^in +^2 K i 



wA.L du 



a(W 2 (y)) a(WHy in )) 
n a(WHy)) n a(Wi(y in )) 



(p»,w») \ , 

du + K^)- 2(r,'e' + ne) - -3(Pi, Wi) 

. '(e 2 ,0) J Z 

+ K 4 (v - v ) + Ksfi(v - v ) , 



(138) 



where W t = y/7ffi) and W^{y) = J^du± J^ i} du — K 



5.4 Solution of the ^-equation 

With (glj} equation ([T2l5j) yields 

?2 n(±l, y)F 2 (±l, y) (±1 - y)H(±l, y)H(y, ±1] 



dt = RE 



(±l-y)H(y,±l) 
This can be rewritten as 



dt 



(l + v±\YG{y) 
Pr(y) dy 



cf^(\I/ + cfl^E) 



dy 



(139) 



(140) 



where Pd{y) are polynomials of order d. P§{y) has the zeros pi (i = 1..6). So we have to solve the integral 

rv Pr{y') dy' 



Mv 1 ) y/W)' 

We apply a partial fractions decomposition where the constants Mi arise: 



t-U 



y^ + M 7 + M 8 y') dy ' 



(141) 



(142) 



If D(±l,y) = 1 (i.e. s — 0) or m = then Y{y) is a polynomial of sixth order and the equation can be 
solved. The hyperelliptic integrals of the first and third kind can be solved as shown in section [3T4l or 14.31 
The solution of the i-equation (|125[) is 



2 



du 



(Pi,Wi, 



du + K 



ln a(W\y)) ln *(WHy in )) 



'(e 2 ,0) J 

+ M 7 (v - vo) + M 8 h{v - vq) 



2(i ? V + t|e)--3(Pi,W i ; 



(143) 



where W t = ^Y{^) and W^ 2 {y) = J^du± J™'" i} du - K x . 



6 The orbits 

One can think of ring coordinates as two pairs of polar coordinates 



x 1 =r lS m(<fi) and x 3 = r 2 sHfi (144) 

X2 = T\ COS(0) X4 = T2 cos(V') 



where 



r 1= R Vl ^ and r 2 =R ^ y2 1 (145) 
x — y x — y 

(see jT5J ) ■ xi, £2, £3, ^4 are four-dimensional Cartesian coordinates. 

If ip is constant, the horizon of the black ring consists of two S 2 spheres. If we look at the rotational axis 
where y = — 1, the coordinates x\ and X2 describe the plane between these two spheres, so the horizons 
and the singularity cannot be seen in this plane. 

If <fi is constant, the horizon has S 1 x S 1 topology. So if x = ±1 the coordinates x$ and X4 describe the 
equatorial plane "as seen from above" , so each horizon consists of two circles. The singularity is a circle 
of radius 1 . 

If both angles are constant, the coordinates describe the x-y-plane. Here we can change from ring coor- 
dinates (x,y) to polar coordinates (/?, 9) via the transformation 



x — y 



Then conventional Cartesian coordinates take the form 



a = p sin 9 , b = p cos 9 



(147) 



(see [H]). The singularity of the black ring is at a = ±1, 6 = 0. 
Note that a = r 2 and b — r\ so that 

x\ — pcos(9) sin(</>) 

x 2 = pcos(6>)cos(» , 14g , 
£3 = psin(9) sm(ip) 
X4 = psin(9) cos(i/>) . 

Figures [TU1 - 1171 show examples for possible orbits in the (charged) doubly spinning black ring spacetime. 
The orbits are shown for cj) = ^ = const. (xi-rz^-rz^-plot) and i\) = 5 = const. (a^-rz^-^s-plot). Sometimes 
the orbit is also plotted in the equatorial plane (x3-a;4-plot) and the plane of constant angles (a-6-plot), 
in these plots we have also marked the ergosphere. 

From the plots where one or both angles were set constant one might get the impression that the 
motion takes place only on one "side" of the black ring. In the x\ — X3 — X4-plots the orbit seems to be 
either above or below the equatorial plane, while in the x\ — X2 — X3-plots the motion seems to take place 
only in the vicinity of one of the two spheres that show the horizon. However, this is an artifact of the 
chosen representation. When none of the angles is set constant, particles and light move all around the 
black ring with respect to their turning points. Nevertheless, we chose to exhibit the orbits with constant 
angle for reasons of clarity. 

A MBO and a TO for light in the ergosphere are depicted in figure \TU\ and [TIJ Figure [12] shows a 
BO for light on the axis of V'-rotation, this orbit also lies in the ergosphere. An EO for particles in the 
equatorial plane of an uncharged black ring can be seen in figure [13J 

A TEO for light around a charged black ring is shown in figure [14J The orbit forms a little loop, which 
can be explained by the fact that there is an ergosphere- free region inside the black ring (see section [2]). 
After the light has left the ergosphere it is no longer dragged along by the rotation of the black ring and 
changes its direction. The light reaches its turning point and approaches again the ergosphere where it 
turns around a second time. 

A MBO for particles in the equatorial plane of an uncharged black ring is depicted in figure [T5] TOs 
in the equatorial plane are shown in figure [16] (particles in the uncharged black ring spacetime) and figure 
IT7I (light in the charged black ring spacetime) . In figure [16] the effect of the ergosphere can be seen as the 
particle changes its direction after entering the ergosphere-free region inside the black ring 

7 Conclusion 

In this paper we presented the analytical solutions of the geodesic equations of the charged doubly spinning 
black ring for special cases. Since the Hamilton- Jacobi equation seems not to be separable in general, we 
had to concentrate on the zero energy null geodesies in the ergosphere (E = m = 0), geodesies on the axis 
of ^-rotation (y = —1) and geodesies on the equatorial plane (x = ±1), which is the axis of ^-rotation. 
We discussed the general structure of the orbits and gave a complete classification of their types. 

In the case E — m — terminating orbits, many- world bound orbits and bound orbits for y = const. = 
— 1 are possible. The charge of the black ring has no effect on the orbits in the ergosphere. If the black 
ring is singly spinning {y = 0) only terminating orbits exist (see [20] ) . 

On the axis of i/)-rotation y is constant, so here the 2;- motion determines the type of orbit. We found 
escape orbits and bound orbits in the (charged) doubly spinning spacetime like in the singly spinning 
black ring spacetime. 

The axis of ^-rotation (the equatorial plane) is divided into two parts: the plane enclosed by the ring 
(x = +1) and the plane surrounding the ring (x — —1). On the plane enclosed by the ring (x = +1) 
terminating orbits and many- world bound orbits are possible. If the doubly spinning black ring is charged, 
bound orbits behind the inner horizon exist. If the black ring is singly spinning only terminating orbits 
are possible. 

On the plane surrounding the ring (x — —1) escape orbits, terminating orbits, two-world escape orbits 
and many-world bound orbits are possible. If the black ring is singly spinning, only escape orbits and 
terminating orbits exist. 

The separability of the Hamilton- Jacobi equation is a coordinate related phenomenon, so one might 
think of a coordinate system in which it would be possible to separate the Hamilton- Jacobi equation in 
general. But recently Igata, Ishihara and Takamori found evidence of chaotic motion in the singly spinning 
black ring spacetime using the Poincare map [30 . From that one could conclude that it is in general not 
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Figure 10: fc = 1, u = 0.1, A = 0.7, $ = 0.4 and * = 0.8 

Many- World Bound Orbit for light in the Ergosphcrc (E = m = 0). The tori (a) or spheres (b) are the inner and outer 
horizons of the black ring. In (c) and (d) the ergosphere is depicted as a light red area with a red dotted border. The 
black and grey dashed circles in (c) and (d) are the inner and outer horizons. 




Figure 11: fc = 1, v = 0.1, A = 0.7, * = -0.2 and * = 3.6 

Termninating Orbit for light in the Ergosphcrc (E — m — 0). The tori (a) or spheres (b) are the inner and outer horizons 
of the black ring. 




Figure 13: v = 0.2, A = 0.9,W = 5, m = 1, s = and E = 1.38 

Escape Orbit for particles in the outer equatorial plane (x — —1 and $ — 0) of an uncharged black ring. The tori (a) or 
spheres (b) arc the inner and outer horizons of the black ring. 
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(c) X3-X4-plot (d) a-6-plot [<p — ip — ^-) 

Figure 14: v = 0.1, A = 0.7,* = 5, m = 0, s = 0.5 and £ = 2 

Two- World Escape Orbit for light in the outer equatorial plane (x — — 1 and $ — 0) of a charged black ring, i.e. x — — 1 
and <I> — 0. The tori (a) or spheres (b) are the inner and outer horizons of the black ring. In (c) and (d) the ergosphcre 
is depicted as a light red area with a red dotted border. The black and grey dashed circles in (c) and (d) are the inner 
and outer horizons. 




Figure 15: v = 0.2, A = 0.9,* = 5, m = 1, s = and E = 0.1 

Many- World Bound Orbit for particles in the outer equatorial plane (cc = — 1 and $ — 0) of an uncharged black ring. 
The tori (a) or spheres (b) are the inner and outer horizons of the black ring. In (c) and (d) the crgosphcrc is depicted 
as a light red area with a red dotted border. The black and grey dashed circles in (c) and (d) are the inner and outer 
horizons. 



(a) xi-Ss-a^-plot (0— ■§■) 



(b) xi-X2-a;3-plot {ip — ^) 



















\ ' /Z''''^ 
— \ — . — 1, / U I : 




\ i- » ^ 


■J I : | I ' ■ 


\ V \ v> ~ 


" ./ /// x 

•TU^ /// 3 


\ \ \ 


/// 






-4- 
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Figure 16: i/ = 0.1, A = 0.8,* = 10, m = 1, s = and E = 3.34 

Terminating Orbit for particles in the outer equatorial plane (x — — 1 and ^> — 0) of an uncharged black ring. The tori 
(a) or spheres (b) are the inner and outer horizons of the black ring. In (c) and (d) the crgosphcrc is depicted as a light 
red area with a red dotted border. The black and grey dashed circles in (c) and (d) are the inner and outer horizons. 
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Figure 17: v = 0.1, A = 0.7,* = 5, m = 0, s = 0.5 and E = 0.1 

Terminating Orbit for light in the inner equatorial plane (x — +1 and $ 
spheres (b) are the inner and outer horizons of the black ring. 



— 0) of a charged black ring. The tori (a) or 



possible to separate the Hamilton- Jacobi equation of singly spinning black rings in any coordinate system. 
It would be interesting to see if chaotic motion also appears in the (charged) doubly spinning black ring 
spacetime. 

The methods shown in this paper can be applied to other black ring spacetimes like for example the 
supersymmetric black ring jS], [5]. This will be done in future work. 
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